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Based on the ultrarelativistic quantum molecular dynamics (UrQMD) model, the effects of the 



density-dependent symmetry potential for baryons and of the Coulomb potential for produced 
mesons are investigated for neutron-rich heavy ion collisions at intermediate energies. The calcu- 
lated results of the A~/A ++ and ir~ /tt + production ratios show a clear beam-energy dependence 
on the density-dependent symmetry potential, which is stronger for the tt~ /ir + ratio close to the 
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pion production threshold. The Coulomb potential of the mesons changes the transverse momen- 
tum distribution of the tt~ /tt + ratio significantly, though it alters only slightly the ir~ and tt + total 
yields. The tt~ yields, especially at midrapidity or at low transverse momenta and the ir~ /tt + ratios 
at low transverse momenta, are shown to be sensitive probes of the density-dependent symmetry 
potential in dense nuclear matter. The effect of the density-dependent symmetry potential on the 
production of both, K° and K + mesons, is also investigated. 
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I. INTRODUCTION 



The study of isospin effects in nuclear matter or finite nuclei is a well established topic 
in nuclear physics, as well as in nuclear astrophysics On the mean field level, this 

means the contribution of the symmetry energy of particles to the total energy (see, e.g., 
the recent Refs. ), and on the two-body collision level, detailed 

experimental data on the free proton-proton and neutron-proton scattering cross sections 
have been obtained as a function of energy ranging from the Coulomb barrier to ultra 
high energy interactions. At low energies (up to several hundreds of MeV/nucleon), large 
differences are visible between the proton-proton and the neutron-proton cross sections. In 
heavy ion collisions (HICs), it is also necessary to study the effect of density on the isospin- 
dependent nucleon-nucleon cross sections. 

Recently, we have explored the density- and temperature-dependence of the nucleon- 
nucleon elastic scattering cross sections in collisions between neutron-rich nuclei at interme- 
diate energies Q. We found a density dependence of the e. a ,,e scattering cross sections. 
This work is based on the theory of quantum hydrodynamics (QHD), in which the interac- 
tion between nucleons is described by the exchange of a, u, p, and S mesons. Since, in the 
QHD theory, both the mean field and the collision term originate from the same Lagrangian 
density, the density-dependence of the symmetry energy is our next subject to investigate. 
Currently the exploration of the isospin dependence of the nuclear interaction ( especially in 
exotic and very isospin-asymmetric systems) has gained high interest 

In particular, the density dependence of the symmetry energy in dense nuclear matter has to 
be studied in detail , be cause it suffers from large uncertainties in the predictions of various 
theoretical models |l5|,ll6|. Based on the QHD theory, with a similar effective Lagrangian 
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18{ |. the main uncertainty of the density dependence of the symmetry en- 



ergy results from the uncertainties of the density- dependent coupling strengths of various 
meson-baryon interactions, especially that of p — N and 5 — N interactions 0. Also, in 
forthcoming experiments at the Rare Isotope Accelerator (RIA) laboratory (USA) and at 
the new international accelerator facility FAIR (Facility for Antiproton and Ion Research) 
at the Gesellschaft fur Schwerionenforschung (GSI, Germany), more neutron-rich beams are 
planned to be adopted, and hence call for further theoretical studies on isospin-dependent 
HICs in the intermediate energy region. 
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In order to obtain detailed information on the density dependence of the symmetry en- 
ergy in dense nuclear matter, quite a few sensitive probes were brought up for experiments 
in recent years: e.g., the tt~/tt + yield ratio, the pion flow, the transverse momentum distri- 
bution of 7T _ /7r + ratios, the neutron-proton differential flow ^, and the isospin equilibration 
and stopping ?| IS, 20]. More recently, we have also investigated l| the sensitivity of the 
ratio on the density-dependent symmetry potentials. Further investigations seem 
to be worthwhile to shed light on the rather unclear isospin-independent and -dependent 
mean-field potentials of the hyperons. Since the contribution of the symmetry energy to the 
whole dynamics of HICs is rather small at Schwerlonen Synchrotron (SIS) energies, it is all 
the more difficult to explore the density dependence of the symmetry energy based on the 
current experimental situation. Hence, in theory, more sensitive probes should be searched 
for. 

Also, such a detailed analysis is quite necessary and important because of the following 
arising questions: Is the effect of the density- dependent symmetry energy on the tt~ /tt + ratio 
influenced by the beam energy? In which momentum, rapidity, or phase-space, region are 
the pion yields more suitable to probe the density- dependent symmetry energy? And, how 
does the Coulomb potential of mesons affect the emission of pions? Gaitanos et al. jf], [?J 
found that for beam energies higher than about 2A GeV, the sensitivity of the tt~ /tt + ratio 
to the form of the symmetry energy at high densities is strongly reduced. Furthermore, they 
found that the ir~ /n + ratio at high transverse momenta pt could offer sound information on 
the density dependence of the symmetry energy. B.-A. Li et al. ^| revisited the n~ /n + ratio 
by using an up-dated, new version of the isospin-dependent Boltzmann-Uehling-Uhlenbeck 
(IBUU04) model, in which an isospin- and momentum-dependent single nucleon potential 
was adopted. They found that the sensitivity of the ir~ /n + ratio on the density-dependence 
of the symmetry energy becomes obvious after considering the momentum-dependent single 
nucleon potential. They also noticed the effect of the Coulomb potential on the 7r production, 
especially on the transverse momentum spectra. Their work implied that it is sufficient to 
measure accurately the low-energy (or the low-transverse-momentum p t ) pions, instead of 
the whole spectrum. In this paper, we address these aspects of the problem on the basis of 
the ultrarelativistic quantum molecular dynamics (UrQMD) model. 

In this work, we also investigate the K°/K + production ratio as a probe of the density- 
dependent symmetry potential in dense nuclear matter. Based on the constituent quark 
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picture the K° meson contains one <i-quark while the K + contains one w-quark, as well as 
one s-quark in each meson. This is closely related to the isospin asymmetry of the nuclear 
medium at SIS energies. Furthermore, experimentally, the effects of system parameters, such 
as beam energy or impact parameter, on the sensitive probes also have to be kept in mind. 
Theoretically, the uncertainty of the isospin- independent nuclear equation of state (EoS), 
as well as the contribution of the Coulomb potential of mesons, might alter the calculated 
results of these probes, hence, they could even modify the conclusions and should be checked 
carefully. 

In the UrQMD model (version 1.3) |3, H 22, Q, adopted here, we find that most 
of the calculations can simultaneously reproduce many experimental measurements, which 
offers a good starting point for studying the isospin effects at SIS energies. In this work, a 
'hard' and a 'soft' Skyrme-type EoS, without momentum dependence, are adopted for central 
(6 = — 2 fm) or near-peripheral (6 = 5 — 8 fm) 132 Sn + 132 Sn reactions at beam energies 
Eb = 0.5A and 1.5A GeV. We consider phenomenological density- dependent symmetry 
potentials, which will be discussed in more detail in the next section. The Coulomb potentials 
of the mesons are switched on or off in order to analyze more conveniently the contributions 
of the various symmetry potentials to the dynamical evolution of the hadrons in the nuclear 
medium. 

The paper is arranged as follows: In Section II, we clarify the inclusion of the isospin- 
dependent part of the mean field in the UrQMD transport model. In Section III, some 
basic isospin effects on the dynamics of the baryons (here the nucleons and A (1232)) in 
SIS energy HICs are shown. In Section IV, firstly the effects of the isospin-independent 
EoS, the beam energy, and the impact parameter on the emitted pion yields and the 7r~ /tt + 
production ratios are discussed. Then, the various phase-space distributions of pion yields 
and the tt~/tt + ratios, with the different density dependences of the symmetry potentials, 
are investigated. At the end of this section, the sensitivity of the K°/K + ratio to the density 
dependence of the symmetry potential is discussed. Finally, the conclusions and outlook are 
given in Section V. 
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II. THE TREATMENT OF THE POTENTIAL UPDATE OF HADRONS IN THE 
URQMD MODEL 



The initialization of neutrons and protons, the corresponding Pauli blocking, the Coulomb 
potential of baryons, and the isospin dependence of the nucleon-nucleon cross sections have 
been introduced explicitly in the standard UrQMD model. In order to study the isospin ef- 
fects in intermediate energy HICs, it is required to further introduce the symmetry potential 
of baryons. 

The isospin-dependent EoS for asymmetric nuclear matter can be expressed as (see, e.g., 



2i) 



e(u, 5) = 6U,S = e Q (u) + e sym (u)5 2 , (1) 

where u = p/po is the reduced nuclear density and 5 = (p n — p p )/p is the isospin asymmetry 
in terms of neutron (p n ) and proton (p p ) densities. The e (u) term is the isospin- independent 
part, which includes the Skyrme and the Yukawa potentials in the UrQMD model 20|. The 
Yukawa parameter is related to the Skyrme parameters since in infinite nuclear matter the 
contribution of the Yukawa potential to the total energy acts like a two-body Skyrme contri- 
bution. For comparison, both a soft (K = 200 MeV) and a hard (K = 300 MeV, the default 
parametrization in this work) EoS are adopted in this work. e sym is the symmetry energy 
per nucleon, in which the kinetic (u^J an d potential (u^J contributions are included. The 
average symmetry potential energy can be expressed as 

eft = e a F(u), (2) 

where e a is the symmetry-potential strength and F(u) is the density-dependent part. If the 
Fermi-gas model is adopted, the symmetry-potential strength e a is related to the symmetry 
energy at normal density, 5*0, by 

S Q ~e a + ^-. (3) 

Here, is the Fermi kinetic energy at normal nuclear density. 

The relativistic mean field calculations p show a small variation of the symmetry kinetic 
energy extracted from various models, and the uncertainty of the symmetry energy with nu- 
clear density results mainly from the symmetry potential energy. Concerning the symmetry 
potential energy, both the symmetry energy coefficient So and the density dependence of the 
symmetry energy are quite uncertain. For So, deduced from the isovector GDR in 208 Pb and 



from the available data of differences between neutron and proton radii for Pb and several 
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Sn isotopes, a rather small range of values is 32 MeV < S < 36 MeV 25 . More recently, 
a new value of S ~ 31 MeV was obtained, based on the consistent folding analysis of the 
p( 6 He, 6 He)p elastic scattering and p( 6 He, 6 Li*)n charge exchange reaction data measured 
at -£/Lab — 41.6 MeV j^j], while Sq = 34 MeV was obtained in a calculation made within 
the relativistic Brueckner framework, using the Bonn A potential (27 1. Thus, it is necessary 
to further investigate the uncertainty of the So value |28| . In this paper, however, the value 
S = 34 MeV is adopted, since we endeavor to investigate the density dependence of the 
symmetry potential energy. 

In order to mimic the strong density dependence of the symmetry potential at high 
densities, we adopt the form of F(u), used in [4], as 

\ F 1 =u l 7 > 
F{u) = \ . (4) 

{F 2 =u-^ a>l 
Here, 7 is the strength of the density dependence of the symmetry potential. We choose 
7 = 0.5 and 1.5, denoted as the symmetry potentials F05 and F15, respectively, a (in F 2 ) 
is the reduced critical density; for u > a, the symmetry potential energy is negative. We 
adopt a = 3 in this paper and the respective symmetry potential is named as Fa3. The 
symmetry potentials F05, F15, and Fa3 are shown in Fig. Q]as a function of the reduced 
nuclear density u, compared with a linear density- dependent symmetry potential. 

From Fig. Q it is apparent that for u < 1, F05>Fa3>F15, whereas for u > 1, 
F15>F05>Fa3. F05 is larger than Fa3 for all densities. This means that for neutron-rich 
HICs at intermediate energies, less neutrons are pushed into the low density region (u < 1) 
as well as into the high density region (u > 1) for density dependence Fa3, as compared to 
F05. In other words, the symmetry-potential parametrization F05 is always stiffer than Fa3 
(at both subnormal and higher densities), except at normal nuclear density. 

Besides nucleons (N), the resonances N*(1440), A(1232), and the hyperons A and S 
should be considered for SIS energy HICs, where the A(1232) is dominating, denoted as A 
in short. For simplicity, the isospin-independent part of the EoS, eo{u), for all other baryons 
is taken to be the same as for the nucleons. The symmetry potentials for the resonances 
are obtained through the constants of isospin coupling (the Clebsch-Gordan coefficients) in 
the process of A [or iV*(1440)] «-> ttN. For hyperons, based on the analysis of the Lane 
potential I29I . I30I ]. we simply take the symmetry potential of hyperons to be nucleon-like. 




FIG. 1: The symmetry potential energy of nucleons as a function of the reduced nuclear density 
u, using various density dependences (see text). 

However, the symmetry potential of excited states of hyperons is not considered, for lack of 
information. 

Combining both, resonances and hyperons, we express the symmetry potential in a unified 
form, which reads as 

^sym ® ^sym & "synn (^) 

where the values of a and (3 for different baryons (B) are listed in Table UJ From this table 
we can see that the symmetry potentials of A~ and S~ are neutron-like and those of A ++ 
and S + are proton-like. On the other hand, the symmetry potentials of A , A + , iV*(1440), 
E°, and A are a mixture of the neutron and proton symmetry potentials. 

In the standard UrQMD model (version 1.3), the cascade mode is usually adopted for 
produced mesons. In order to investigate the phase-space distributions of pions, especially 
the transverse momentum distributions in intermediate energy HICs, the Coulomb potential 
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TABLE I: Values a and [5 for the symmetry potentials (Eq.|3J) of different nucleon and A resonances, 
as well as hyperons. 
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N*°(1440) 


1/3 2/3 
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1/2 1/2 


N*+(1440) 


2/3 1/3 
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E~ 
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2/3 1/3 




1/2 1/2 


A+ 


1/3 2/3 


s+ 


1 


A++ 


1 



of mesons should be considered explicitly |fj, • Thus, in this paper, we also investigate 
the contribution of the Coulomb potentials between mesons and baryons. Other mean-field 
potentials of mesons are not yet considered. 



III. EFFECTS OF THE SYMMETRY POTENTIAL ON THE DYNAMICS OF 
NUCLEONS AND A'S 

First of all, let us discuss the basic consequences of the density dependence of the sym- 
metry potentials in neutron-rich HICs at SIS energies. Fig. El (top) shows the time evolution 
of the neutron and proton numbers for the symmetry potentials F15, F05, and Fa3. Here, 
we do not distinguish whether the nucleons are free or bound in heavier fragments. The 
beam energies are = 0.5A and 1.5A GeV and central collisions (6 = — 2 fm) are chosen. 
From this plot, we notice that independent of the choice of the density dependence of the 
symmetry potential, the neutron numbers reach their minima at ~ 15 fm/c and 10 fm/c, for 
the beam energies = 0.5A and 1.5A GeV, respectively. At this time the central nuclear 
density reaches its maximum and some of the neutrons are excited (to N* or A) or converted 
to other baryons (e.g., hyperons) through various baryon-baryon or meson-baryon collisions. 
We also notice that the higher the beam energy, the more excited baryon states are pro- 
duced. There are obviously less excited protons than neutrons at the compression stage. In 
other words, at the compression stage, the effect on the neutron number is stronger than on 
the proton number, which is due to the neutron-rich nuclear environment. After the decay of 
these unstable baryons, the emitted mesons (especially, the pions) alter the neutron/proton 
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ratio (the neutron number decreases, and, correspondingly, the proton number increases), 
and hence the ratio tends to become unity with an increase of the beam energy. In addition 
the effect of the density-dependent symmetry potentials on the dynamical process is smaller 
at higher energies JjJ. With a softer symmetry potential Fa3, more neutrons are excited 
at the compression stage, which is clearly due to more neutrons being kept by the softer 
symmetry potential in the high density region, as implied from Fig.^ 

In Fig. |21 (bottom), the neutron/proton ratio of the all (free and bound) nucleons at freeze- 
out time (tf = 50 fm/ c, the maximum time used in the top plot) is shown as a function of the 
normalized rapidity (y^ = y c /ybca,m, where y c is the rapidity of the particle in the center-of- 
mass system and yb e &m is the beam rapidity). At E h = 0.5A GeV, the density dependence of 
the symmetry potential influences the neutron/proton ratio strongly, whereas at = 1.5 A 
GeV the sensitivity is reduced, especially at midrapidity. This happens not only due to the 
relatively weak effect of the EoS, as compared to the two-body collision dynamics, but also 
due to the nucleon-nucleon cross sections at higher energies that depend less on isospin. 
Secondly, in the projectile and target rapidity regions, the n/p ratio is more sensitive to the 
density dependence of the symmetry potential as compared to the midrapidity region. Less 
collisions take place in the projectile and target rapidity regions than in the midrapidity 
region (see also, Fig. |3J). Thirdly, the n/p ratio is larger with a softer symmetry potential 
Fa3 than that with F15, and vice versa in the projectile-target regions. The nucleons in the 
midrapidity region represent mainly the behavior of symmetry potential at high densities, 
where more neutrons are kept with a soft symmetry potential. We should also notice that at 
midrapidity, the difference of the ratios with the F05 and Fa3 symmetry potentials is almost 
negligible. From Fig. E]we can see that F05 is always stiffer than Fa3 at both low and high 
densities, thus the n/p ratio is also influenced by the density dependence of the symmetry 
potential at subnormal densities. Therefore, in order to investigate any isospin dependences, 
it is very important to know the correct density- dependent form of the symmetry potential 
at both, low and high densities. 

Fig. El shows the calculated rapidity distribution of the collision number of neutrons and 
protons (upper plot), and the corresponding ratios (lower plot). Evidently, at both energies, 
the collision number of nucleons has a strong maximum at midrapidity. The collision number 
increases with beam energy, especially for neutrons because of the neutron-rich system. The 
average proton collision number is always larger than the average neutron collision number 
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FIG. 2: Top: Time evolution of the neutron and proton numbers for central Sn+ Sn collisions 
at beam energies 0.5^4 and 1.5^4 GeV. For the "hard" EoS, three different forms of the symmetry 
potential are used (see text). The two horizontal solid (thick) lines represent the initial neutron and 
proton numbers of the system. Bottom: normalized rapidity distribution of the neutron/proton 
ratio of all nucleons at freeze-out time (see text). The error bars show the statistical error. 
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FIG. 3: Top: Normalized rapidity distributions of the collision number of neutrons and protons 
at energies E^ = 0.5A and 1.5A GeV. The case of the symmetry potential F15 is shown. Bottom: 
The corresponding ratios between the collision numbers of neutrons and protons. 



11 



because of the differences in the nucleon-nucleon cross sections at intermediate energies 



11 1 . When the beam energy increases from 0.5A to 1.5A GeV, this isospin effect of the 
nucleon-nucleon cross section is largely reduced, so that the ratio between the neutron and 
proton collision numbers approach unity, which is shown in the lower part of Fig. El In the 
midrapidity region, the collision ratio of neutrons and protons is ~ 82% at = 0.5A GeV, 
which increases to ~ 95% at E\, = 1.5 A GeV. 

Since, in the SIS energy region, 7r-mesons are mainly produced from decaying A's, it is 
also important to investigate the dynamics of the A's with respect to the symmetry potential. 
Fig. |3] shows the time evolution of the various components of A (upper plot) and the ratio 
A~/A ++ (lower plot) for the symmetry potentials F15 and Fa3 and for central 132 Sn+ 132 Sn 
collisions at the two energies E\> = 0.5 A and 1.5 A GeV. In the upper plot, only the case 
£b = 0.5A GeV is shown since, at E^ = 1.5A GeV, the effect of the density dependence of 
the symmetry potentials on the A production is almost negligible. This is supported by the 
time evolution of the ratio A~/A ++ at E\> = 1.5A GeV in the lower plot. 

From the upper part of Fig. 01 it is clear that the effect of the density-dependent symmetry 
potentials on the time evolution is strongest for the A~ production, like for neutrons in 
Fig. El From the time evolution of the A~/A ++ ratio in the lower plot of Fig. 01 we further 
notice that, similar to the case of nucleons, the effect of the density-dependent symmetry 
potentials on the A~/A ++ ratio is largely reduced at the higher beam energy 1.5A GeV. 
And the time evolution of this ratio is nearly flat since, at this energy, the effect of the 
Coulomb and the symmetry potentials becomes small. At the lower beam energy 0.5A GeV 
and during the compression stage (t < 12 fm/c), the A~/A ++ ratio decreases with time 
while it increases at the later stage, the expansion stage. This is apparently due to the 
dynamically mutual interaction between the Coulomb and the symmetry potentials: along 
the time evolution at the compression stage, more neutrons are pushed out because the 
relative strength of the symmetry potential as compared to the Coulomb potential increases. 
This can be understood from the study on the ratio of preequilibrium neutron number to 
proton number in the intermediate-energy neutron- rich HICs in Ref. j^ ]: this ratio is larger 
than the initial neutron/proton ratio of the colliding system; at the expansion stage, more 
protons are pushed out due to the stronger Coulomb potential as compared to the symmetry 
potential. Note that the increase of the A~/A ++ ratio at the late stage is stronger for the 
symmetry potential Fa3 than for F15. The time evolution of the A abundancies, depending 
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on the symmetry potential, should also affect the production of pions. 

Besides the 7r-yields, the pion spectra should also be influenced by the decay of A's. Fig.|3] 
shows the transverse momentum spectrum of A's (upper plot) and the A~/A ++ ratios (lower 
plot) at t = 20 fm/c and = 0.5A GeV, for the symmetry potentials F15, F05, and Fa3. 
We notice from the upper part of Fig. EJthat, similar to the neutrons in Fig. |2J the effect 
of different density-dependent symmetry potentials is strongest for the A - . From the lower 
plot of Fig. El we can see that, at lower p£ m (pj m < 0.8 GeV/c), the A~/A ++ ratios are 
quite different for the symmetry potentials F15 and F05 (or Fa3), and the results of the 
F05 and Fa3 symmetry potentials are nearly the same. These results are similar to the 
rapidity distribution of the neutron/proton ratio, shown in Fig- El A's with low transverse 
momenta are not produced from very high densities, namely around normal density, where 
the difference between F05 and Fa3 is quite small and the effect of F05 and Fa3 on the 
A~/A ++ ratio are cancelled strongly between low and high densities. For p£ m > 0.8 GeV/c, 
the differences between the results of F05 and Fa3 become distinguishable because A's with 
the high transverse momenta are mainly produced from the central high-density region. 
Thus the A~/A ++ ratio at high transverse momenta can provide clearer information on 
the density dependence of the symmetry potential. In the next section, we investigate 
the question to what extent this property could be transferred to tt's, especially when the 
contribution of the Coulomb potential of mesons is also considered. 

IV. EFFECTS OF THE SYMMETRY POTENTIAL ON THE vr AND K PRODUC- 
TION 

Fig. IH1 shows the time evolution of the 7r abundancies (upper plot) and their ratios (lower 
plot) for the density-dependent symmetry potentials F15, F05, and Fa3. Here, the Coulomb 
potential of mesons is switched off. However, for comparison, we also provide the results 
(for the case of F05 and at = 0.5A GeV) when the Coulomb potential of mesons is 
taken into account (solid points in Fig. El at time t = 45 fm/c). Apparently, the effect 
of the Coulomb potential of mesons on the 7r-yields is quite small. We also find that the 
effect of the different symmetry potentials on the tt~/tt + ratio is hardly affected by the 
Coulomb potential of mesons although it alters the absolute value of the ratio. The number 
of 7r~ grows much stronger than the number of 7r + , almost independent of the choice of the 
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symmetry potentials F05 and Fa3. The 7r~-yields are, however, smaller for F15. So is the 
7r~/7r + ratio at = 0.5A GeV. On the other hand, for the higher beam energy E\> = 1.5A 
GeV, the ratios are almost independent of the symmetry potentials. (These results are 
clearly due to the properties of A's, shown in Figs. H]and[SJ). Therefore, it is necessary to 
pay attention to the energy dependence of the isospin effect on the ir" /n + ratios as well. 

Next, we show the dependence of the tt~/tt + ratio on various isospin-independent EoS 
and impact parameter b. In view of the above noted beam-energy dependence, we only 
consider the lower beam energy = 0.5A GeV in the following. Fig. [7| shows the results 
for two EoS and two centralities. The effect of the density-dependent symmetry potentials 
on the tt~/tt + ratios is almost not affected by the uncertainty of the isospin-independent 
EoS. That is, the absolute shift in the n~/7r + ratio due to different symmetry potentials 
remains almost the same, independent of the stiffness of the isospin-independent EoS. On 
the other hand, at larger impact parameters (b = 5 ~ 8 fm), the tt~/tt + ratio is even more 
sensitive to the density dependence of the symmetry potential. 

The rapidity y c and the polar angle 9 C (in the center-of-mass system) distributions of 
7r~ and 7r + for various density-dependent symmetry potentials are shown in Fig. El The 
influence of the density-dependent symmetry potentials on the 7r-multiplicity distributions 
is strongest at midrapidity and at polar angles around 90°. This holds in particular for the 
7r~ distribution, which was also seen in Q for the kinetic energy spectrum of pions in central 
132 Sn+ 124 Sn reactions at E h = OAA GeV. 

The upper two plots of Fig. El show the transverse momentum distributions of tt~ and 
7r + , with and without the contribution of the Coulomb potentials of mesons, as well as with 
a rapidity-cut (\y c \ < 0.2) (left plots) and without any kinematical cut (right plots), for the 
symmetry potential F05. This is important because, in recent years, several experiments on 
charged pion production at intermediate energies were performed by the FOPI Collaboration 
at GSI where the possible contribution of the Coulomb potentials of pions to the 7r-transverse 
momentum distribution was also discussed j^jH^. Although, as shown in Fig. [HI the effect 
of the Coulomb potential of mesons on the total ir yield is small, Fig. El demonstrates that 
it strongly influences the momentum distribution; the Coulomb potential of mesons always 
shifts the 7r~ to lower, and the ir + to higher p£ m for the neutron-rich systems, in line with 
Ref. 0|. With the rapidity-cut \y c \ < 0.2, this effect becomes even more pronounced; as a 
result, the transverse momentum distribution of the n~ /tt + ratio at low pf 11 becomes steeper 
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than the one without the rapidity cut (lower plots in Fig. El show the corresponding tt~/tt + 
ratios for the symmetry potentials F15, F05, and Fa3). If the Coulomb potential of mesons 
is switched off, the tt~/tt + ratio without rapidity-cut is nearly constant at low transverse 
momenta (p£ m < 0.3 GeV/c). With the rapidity-cut (\y c \ < 0.2) it increases weakly with 
p£ m because of the contribution of the Coulomb potential of A's. At higher transverse 
momenta (p t > 0.4 GeV/c), and with rapidity-cut (lower- left plot of Fig. EJ), the effect of 
the density-dependent symmetry potentials on the n~ /n + ratios is noticeable; the difference 
is also obvious for F05 and Fa3. It is recalled that a similar phenomenon can be seen in 
Fig. El for the A~/A ++ ratios. However, when the Coulomb potentials of mesons are taken 
into account, the effect of the density- dependent symmetry potentials on the 7t~/tt + ratios 
at high transverse momenta is largely reduced, regardless of a rapidity cut. It was pointed 
out in that the Coulomb potential is stronger than the symmetry potential, and the 
Coulomb potential acts directly on the charged pions but the symmetry potential does not, 
such that the 7r's emitted from the dense region, which have higher p£ m , may be influenced 
much more by their Coulomb potential than by the symmetry potential. Considering the 
limit of experiments so far, we claim that the 7r~/7r + production ratio at p^ m = 0.1 GeV/c is 
a suitable candidate for probing the density dependence of the symmetry potential, which 



is discussed in more detail in 



Mi- 



Fig. El shows the density distribution of pion emission, with and without Coulomb po- 
tential of mesons. In the upper plot, the percentages of each component of ir and of the 
total it yield are shown as a function of the reduced density u, in steps of Au = 0.2, with the 
Coulomb potential of mesons included. Most pions are emitted from densities higher than 
normal nuclear density (especially from the density region u ~ 1 — 2); a few pions (~ 16%) 
are emitted from subnormal densities. This follows from Fig. El where quite a few A's are 
transported close to normal nuclear density. As these A's decay, some of them also appear 
at subnormal densities. 

In the lower part of Fig. EH the ratios of the percentages of 7r~ and 7r + are shown as 
a function of the reduced density u, with and without the Coulomb potential of mesons. 
The ratios are always larger at lower densities than at higher densities since at the late 
stage, where the nuclear density decreases, the A~/A ++ ratio (shown in Fig. 0]) increases 
with time. When the Coulomb potential of mesons is taken into account, more ir~ mesons 
are decelerated because of the larger amount of the negatively charged particles in the low 
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density region, and are thus emitted with low momenta which leads to the rise of the ir~ /n + 
ratios at low p£ m , as shown in Fig. El 

Besides the tt~ /h + ratios, the K°/K + ratios are also thought to be a suitable candidate 
to probe the density dependence of the symmetry potential in dense nuclear matter, as 
already discussed in the Introduction. Fig. ^2 shows the time evolution of the K° and 
K + abundancies, as well as their ratios K°/K + , for central 132 Sn+ 132 Sn collisions at beam 
energy 1.5A GeV and for the symmetry potentials F15 and Fa3. Similar to the 7r production, 
the effect of different density dependences of the symmetry potential on the kaon yields is 
quite reduced at higher energies. 

In SIS energy HICs, kaons are emitted mainly from the high density region and at the 
early stage of the reaction. The dominant production channels are BB —>■ BB* — ► BKY 
and Bit — > B* — > KY (here Y represents a hyperon). With a softer symmetry potential, 
in the neutron-rich nuclear medium, more neutrons (protons) are shifted to the high (low) 
density region and hence more negatively (less positively) charged B* are produced, and 
the K°/K + ratio increases. Thinking of the energy dependence, the kaon production at 
energies much lower than its threshold (~ 1.58 GeV from the nucleon-nucleon interaction in 
free space) might be more sensitive to the density dependence of the symmetry potential. 
As an example, we have also calculated the kaon yields from the reaction 208 Pb + 208 Pb at 
£b = 0.8A GeV and b = 7 ~ 9 f m with the symmetry potentials F15 and Fa3: the K°/K + 
ratio for F15 is about 1.25, whereas it is about 1.4 for Fa3. We should point out that these 
results on kaons might be relatively rough since we do not consider any kaon-nucleon mean- 
field potential [3, lsfil 3t| . As far as we know, there is no explicit calculation about the 
difference of K + - and i^°-nucleon potentials in the nuclear medium. 



V. SUMMARY 



In this paper, we have investigated the role of the density-dependent symmetry potential 
in heavy ion collisions (HICs) at SIS energies, based on the UrQMD model (vl.3). The 
contribution of the Coulomb potential of mesons has been studied, in order to better un- 
derstand the effects of the density-dependent symmetry potential on the dynamics of pions 
produced in neutron-rich HICs. The calculated results show a strong beam-energy depen- 
dence of the effect of the density-dependent symmetry potentials on the production ratios 



16 



A~/A ++ as well as n~ /n + . The uncertainty in the isospin-independent EoS alters the values 
of both the pion yields and the 7r~/7r + ratios, but has almost no effect on the role of the 
density-dependent symmetry potential on the n~ /n + ratios. The impact parameter is found 
to be important for the 7r~/ii + ratios; for a larger impact parameter the effect of different 
density-dependent symmetry potentials on the ir~ /n + ratios becomes stronger. 

The Coulomb potential of mesons changes the transverse momentum distribution of the 
ty~/'k + ratios significantly, though it leaves the total n~ and ir + yields almost unchanged. 
We find that the negatively charged pion yields, especially at midrapidity and low transverse 
momenta, as well as the n' /n + ratios at low transverse momenta, could be sensitive to the 
density-dependent symmetry potential in a dense nuclear matter. Quite a few pions are still 
produced at subnormal densities and thus are affected by the density-dependent symmetry 
potential at subnormal densities. These studies are of interest for forthcoming experiments 
at RIA (USA) and FAIR/GSI (Germany). 

In this work, we have also investigated the yields of K° and K + mesons and the ratios 
K°/K + from neutron-rich HICs at E h = 1.5 A GeV. It is shown that they do not seem to be 
suitable to investigate the density dependence of the symmetry potential in dense nuclear 
matter. However, kaon production at energies much lower than threshold might improve 
the situation and is worth investigating further. Meanwhile, it is necessary to consider the 
difference of the K + - and i^°-nucleon potentials in the nuclear medium. 
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FIG. 4: Top: Time evolution of A abundancies, for the symmetry potentials F15 and Fa3 at the 
beam energy = 0.5^4 GeV. Bottom: Time evolution of the A~/A ++ ratio for the two symmetry 
potentials and at two beam energies E^ = 0.5A and 1.5^4 GeV. 
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FIG. 5: Top: Transverse momentum distributions of various A components for the symmetry 
potentials F15, F05, and Fa3, at E\, = 0.5 A GeV, and t = 20 fm/c. Bottom: The corresponding 
A~/A++ ratios. 
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FIG. 6: Time evolution of the tt abundancies (top) and the 7r — /vr + ratios (bottom) for the symmetry 
potentials F15, F05, and Fa3, at the beam energy = 0.5^4 GeV. The Coulomb potential of 
mesons is switched off, except for the solid dots for the case F05 at E\, = 0.5A GeV. In the lower 
plot, the 7r~ /tt + ratio at E^ = 1.5A GeV is also shown. 
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FIG. 7: 7r /ir + ratios from 132 Sn+ 132 Sn reactions at 0.5^4 GeV, using a 'hard' and a 'soft' EoS 
with impact parameters b = — 2 fm and 5 — 8 fm. The symmetry potentials used are F15 and 
Fa3. 
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FIG. 8: Rapidity (top) and the polar angle (bottom) distributions of n and ir + in central 
132 Sn+ 132 Sn collisions at = 0.5A GeV, for the symmetry potentials F15 and Fa3. 
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FIG. 9: Upper two plots: Transverse momentum distribution of ir~ and tt + mesons with (left) and 
without (right) rapidity cut \y c \ < 0.2. The results with and without Coulomb potential of mesons 
(denoted as "M Coul.") are shown. The symmetry potential used is F05. Lower two plots: The 
corresponding tt~ /n + ratios, but for three different density-dependent symmetry potentials. 
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FIG. 10: Top: Percentages of the 7r's as a function of reduced density, with the Coulomb potential 
of mesons included. Bottom: the corresponding density distribution of the 7r~ /tt + ratio of the 
percentages, with or without the Coulomb potential of mesons. 
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FIG. 11: Top: Time evolution of the K° and K + abundancies for central 132 Sn+ 132 Sn collisions at a 
beam energy 1.5^4 GeV and with the symmetry potentials F15 and Fa3. Bottom: the corresponding 
time evolution of the K°/K + ratios. 
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